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Abstract 


We  give  an  O {y/nm  log  N)  algorithm  for  the  single-source  shortest  paths  problem  with 
integral  arc  lengths.  (Here  n  and  m  is  the  number  of  nodes  and  arcs  in  the  input  network 
and  N  is  essentially  the  absolute  value  of  the  most  negative  arc  length.)  This  improves 
previous  bounds  for  the  problem. 


1  Introduction 

In  this  paper  we  study  the  shortest  paths  problem  where  arc  lengths  can  be  both  positive  and 
negative.  This  is  a  fundamental  combinatorial  optimization  problem  that  often  comes  up  in 
applications  and  as  a  subproblem  in  algorithms  for  many  network  problems.  We  assume  that 
the  length  function  is  integral,  as  is  the  case  in  most  applications. 

We  describe  a  framework  for  designing  scaling  algorithms  for  the  shortest  paths  problem  and 
derive  several  algorithms  within  this  framework.  Our  fastest  algorithm  runs  in  0{y/nm\ogN) 
time,  where  n  and  m  are  the  number  of  nodes  and  arcs  of  the  input  network,  respectively,  and 
the  arc  costs  are  at  least  -N}  Our  approach  is  related  to  the  cost-scaling  approach  for  the 
minimum-cost  flow  problem  [2,  11,  13,  14,  15]. 

Previously  known  algorithms  for  the  problem  are  as  follows.  The  classical  Bellman-F^d 
algorithm  [l‘ 7]  runs  in  0{nm)  time.  Our  bound  is  better  than  this  bound  for  N  =  o{2'^). 
The  scaling  algorithm  of  Gabow  and  Tarjan  [10]  is  dominated  by  an  assignment  subroutine, 
which  takes  Oiy/Tivi  log(nC'))  time,  where  [-C, . . . ,  C]  is  the  range  of  the  arc  costs.  Our  bound 
dominates  this  bound  (note  that  <  C  —  1). 

Our  bound  is  competitive  even  with  bounds  for  special-purpose  algorithms  on  planar 
graphs.  The  fastest  shortest  paths  algorithm  currently  known  for  planar  graphs  [8]  runs  in 
0(7i''^logu)  time.  Our  algorithm  runs  in  O(n^'^logiV)  time  on  planar  graphs,  which  is  bettei 

for  N  =  o{n). 

The  previous  best  time  bounds  for  scaling  algorithms  for  the  shortest  paths  problem  match 
those  for  the  assignment  problem.  Our  improved  bound  for  the  former  problem  implies  that 
either  this  problem  is  computationally  simpler  or  the  bound  for  the  latter  problem  can  be 
improved. 

Our  framework  is  very  flexible.  In  Section  8,  for  example,  we  use  the  framework  to  design 
a  different  algorithm  for  the  problem  that  runs  in  log(niV))  time.  Although  our  algo¬ 

rithms  do  not  use  sophisticated  data  structures  such  as  those  for  manipulating  sets  and  paths, 
the  flexibility  of  the  framework  would  make  it  easy  to  use  such  data  structures  if  needed.  This 
flexibility  may  lead  to  better  running  time  i)ounds. 


2  Preliminaries 

The  input  to  the  single-source  shortest  paths  problem  is  (6'.  -s,  /),  where  G  =  (F,  E)  \s  a;[''‘^ected 
graph,  f  :  E  R  is  a  length  function,  and  .s  €  V  is  the  source  node.  (See  e.g.  [3,  !<]•)  ^  “e 

'  We  assume  that.  iV  >  2  so  that  log  A  >  0. 
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goal  is  to  find  shortest  paths  distances  from  s  to  all  other  nodes  of  G  or  to  find  a  negative 
length  cycle  in  G.  If  G  has  a  negative  length  cycle,  we  say  that  \  \e  problem  is  infeasible.  We 
assume  that  the  length  function  is  integral.  We  also  assume,  without  loss  of  generality,  that 
all  nodes  are  reachable  from  s  in  G  and  that  G  has  no  multiple  arcs.  The  latter  assumption 
allows  us  to  refer  to  an  arc  by  its  endpoints  without  ambiguity. 

We  denote  1^1  by  n  and  lEj  by  m.  Let  M  be  the  smallest  arc  length.  Define  N  =  -M  if 
M  <  —  1  and  N  =  2  otherwise.  Note  that  N  >2  and  /(a)  >  -N  for  all  a  £  E, 

A  price  function  is  a  real-valued  function  on  nodes.  Given  a  price  function  p,  we  define  a 
reduced  cost  function  /p  ;  £*  R  by 

lp{v,  lu)  =  /(u,  iy)  +  p{v)  -  p{w). 

We  say  that  a  price  function  p  is  feasible  if 

/p(a)  >0  \/aeE,  (1) 

For  an  €  >  0,  we  say  that  a  price  function  is  e^feasible  if 

lp{a)  >  \/aeE,  (2) 

We  call  an  arc  (u,ty)  improvable  if  /p(y,ty)  <  —€,  and  we  call  a  node  w  improvable  if  there 
is  an  improvable  arc  entering  iv. 

Given  a  price  function  p,  we  say  that  an  arc  a  is  admissible  if  lp{a)  <  0,  and  denote  the  set 
of  admissible  arcs  by  Ep.  The  admissible  graph  is  defined  by  Gp  =  (V,  Ep). 

If  the  length  function  is  nonnegative,  the  shortest  paths  problem  can  be  solved  in  0(m  + 
nlogT?)  time  [9].  We  call  such  a  problem  Dijkstra's  shortest  paths  problem  [5].  Given  a  feasible 
price  function  p,  the  shortest  paths  problem  can  be  solved  as  follows.  Let  d  be  a  solution  to 
the  Dijkstra’s  shortest  paths  problem  (G,5, /p).  Then  the  distance  function  d'  defined  by 
(P[v)  =  d[v)  +  p{v)  -  p(5)  is  the  solution  to  the  input  problem. 

We  restrict  our  attention  to  the  problem  of  computing  a  feasible  price  function  or  finding 
a  negative  length  cycle  in  G. 


3  Successive  Approximation  Framework 


Our  method  computes  a  sequence  of  €-feasible  price  functions  with  (  decreasing  by  a  factor  of 
two  at  each  iteration.  Initially,  all  the  prices  are  zero  and  f  is  the  smallest  power  of  two  that  is 
greater  than  N.  The  method  maintains  integral  prices.  At  each  iteration,  the  method  halves 
e  and  applies  the  refine  subroutine,  which  takes  as  input  a  (2£)-feasible  price  function  and 
returns  an  {-feasible  price  function  or  discovere  a  negative  length  cycle.  In  the  latter  case,  the 
computation  halts.  The  high-level  description  of  the  method  is  given  in  Figure  1. 

Lemma  3.1  Suppose  n  price  function  p  is  integral  and  l-feasible.  Then  for  every  a  £  E. 
Cp(a)  >  0. 
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procedure  C0MPUTE-PR1CES(G,  s, /); 
[initialization] 

e  4_ 

Vt)  6  E,  p{v)  <-  0; 

[loop] 

while  €  >  1  do 
(  e/2; 

p  REFlNE(e,p); 

end; 

return  (p); 

end. 


Figure  1:  High-level  description  of  the  shortest  paths  method. 

Proof.  The  lemma  follows  from  the  fact  that  ip{a)  is  integral  and  ip{a)  >  -1.  ■ 

Corollary  3.2  The  method  terminates  in  0{\ogN)  iterations. 

4  Dealing  with  Admissible  Cycles 

Suppose  that  Gp  has  a  cycle  F.  Since  the  reduced  cost  of  a  cycle  is  equal  to  the  length  of  the 
cycle,  /(P)  <  0. 

If  /(F)  <  0,  or  /(F)  =  0  and  there  is  an  arc  (v,  iu)  such  that  lp{v,  -w)  <  0  and  both  v  and 
(u  are  on  F,  then  the  input  problem  is  infeasible  and  the  method  terminates.  Otherwise,  we 
contract  F  and  remove  self-loops  adjacent  to  the  contracted  node.  .A.  feasible  price  function  on 
the  contracted  graph  extends  to  a  feasible  price  function  on  the  original  graph  in  a  straight- 
forward  way. 

Our  algorithm  uses  an  0(m)-time  subroutine  DECYCLE(GV)  that  works  as  follows.  Find 
strongly  connected  components  of  Gp  [1C];  if  a  component  contains  a  negative  reduced  cost 
arc,  G  has  a  negative  length  cycle;  otherwise,  contract  each  component. 

Suppose  Gp  is  acyclic.  Then  Gp  deHnes  a  partial  order  on  V  and  on  the  subset  of  improvable 
nodes.  This  motivates  the  following  definitions.  A  set  of  nodes  5  is  closed  if  every  node 
reachable  in  Gp  from  a  node  in  S  belongs  to  S.  A  set  of  nodes  (arcs)  5  is  a  chain  if  there  is  a 
path  in  Gp  containing  every  element  of  S. 

5  Cut-Relabel  Operation 

In  this  section  we  study  the  CUT-REL.ABEL  operation  which  is  used  by  our  method  to  transform 
a  (2€)-fea.sible  price  function  into  an  e-feasible  one.  The  cut-relabel  operation  takes  a  closed 
set  5  and  decreases  prices  of  all  nodes  in  5  by  e. 

Lemma  5.1  The  CUT-REL.-\BEL  operation  does  not  create  any  improtHible  arcs. 
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Proof.  The  only  arcs  whose  reduced  cost  is  decreased  by  cut-relabel  are  the  arcs  leaving 
S.  Let  a  be  such  an  arc.  The  relabeling  decreases  Zp(a)  by  e.  Before  the  relabeling,  5  is  closed 
and  therefore  lp{a)  >  0.  After  the  relabeling,  Zp(n)  >  ■ 

The  above  lemma  implies  that  cut-RELABEL  does  not  create  improvable  nodes.  The  next 
lemma  shows  how  to  use  this  operation  to  reduce  the  number  of  improvable  nodes. 

Lemma  5.2  Let  p  be  a  {‘li) -feasible  price  function.  Let  S  be  a  closed  set  of  nodes,  and  let 
X  C  S  be  a  set  of  improvable  nodes  such  that  every  improvable  arc  entering  X  crosses  the  cut 
defined  by  S.  After  the  set  S  is  relabeled,  nodes  in  X  are  no  longer  improvable. 

Proof.  Let  p'  be  the  price  function  after  the  relabeling.  Let  w  e  X  and  let  (v,tu)  be  an 
improvable  arc'with  respect  to  p.  By  the  statement  of  the  lemma,  v  ^  S.  Thus  the  relabeling 
increases  Ip  by  £,  and,  by  (2€)-feasibility  of  p,  lpi{v,x)  >  —e.  ■ 

A  simple  algorithm  based  on  cut-relabel  applies  the  following  procedure  to  every  im¬ 
provable  node  V. 

1.  DECYCLE(Gp). 

2.  5  set  of  nodes  reachable  from  {u}  in  Gp. 

3.  cut-relabel(5). 

4.  If  t;  is  improvable  return  “the  problem  is  infeasible”  and  halt. 

It  is  easy  to  see  that  given  a  (2£)-feasible  price  function,  this  algorithm  computes  an  £-feasible 
one  in  0{nm)  time. 

As  we  shall  see,  it  is  possible  to  find  either  a  set  A',  such  that  relabeling  A'  eliminates  many 
improvable  nodes,  or  a  chain  containing  many  improvable  arcs.  In  the  next  section  we  describe 
a  technique  that  can  be  applied  to  a  chain  of  improvable  arcs. 

6  Eliminate- Chain  Subroutine 

Suppose  that  Gp  is  acyclic  and  let  F  be  a  path  in  Gp.  Let  (vi,Wi), . . . ,  (v/,  wt)  be  the  collection 
of  all  improvable  arcs  on  F  such  that  for  1  <  t  <  j  <  t,  the  path  visits  i?,  before  n,  (t.c..  i  i 
is  visited  last).  By  definition,  nodes  are  improvable.  In  this  section  we  describe 

a  subroutine  ELIMINATE-CHAIN  that  modifies  p  so  that  the  nodes  wi,...,u’t  are  no  longer 
improvable  and  no  new  improvable  nodes  are  created,  or  finds  a  negative  length  cycle  in  G. 
The  subroutine  runs  in  0{m)  time. 

At  iteration  i,  ELIMINATE-CHAIN  finds  the  set  5,  of  all  nodes  reachable  from  tn,  in  the  ad¬ 
missible  graph  and  relabels  S';.  If  tu;  is  improvable  after  the  relabeling,  the  algorithm  concludes 
that  the  problem  is  infeasible. 
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Lemma  6.1  The  path  T  is  always  admissible.  If  Wi  is  improvable  after  iteration  i,  then  the 
problem  is  infeasible. 

Proof  The  price  function  is  modified  oniy  by  cUT-RELABEL.  At  iteration  i,  5i  contains  w,,  ail 
its  successors  on  T,  and  no  other  nodes  of  F  (by  induction  on  i).  Therefore  (,(<,(, »,)  changes 
exactly  once  during  iteration  i,  when  it  increases  by  t.  The  arc  (now)  U  improvable  before 
the  change,  and  admissible  after  the  change.  Reduced  costs  of  other  arcs  on  T  do  not  change 
during  the  execution  of  eliminatE-CHAIN. 

Suppose  is  improvable  immediately  after  iteration  t.  Then  there  must  be  a  node  u  such 
that  (v  Wi)  is  improvable  and  u  G  5^.  By  construction  of  S.,  there  must  be  an  admissible  path 
from  Wi  to  u.  This  path  together  with  the  arc  {v,  Wi)  forms  a  negative  length  cycle.  ■ 

Lemmas  5.1  and  6.1  imply  that  the  implementation  of  ELIMINATE-CHAIN  is  correct.  Next 
we  show  how  to  refine  this  implementation  to  achieve  0(m)  running  time.  The  key  fact  that 
allows  such  an  implementation  is  that  the  sets  Si  are  nested. 

First  we  contract  the  set  of  nodes  at  every  iteration.  The  reason  for  contracting  is  to 
allow  us  to  change  the  prices  of  nodes  in  Si  efficiently  (these  prices  change  by  the  same  amount) 
The  coNTRACT(Si)  operation  collapses  all  nodes  of  Si  into  one  node  s.-  assigns  the  price  of 
the  new  node  to  be  zero,  and  eliminates  self-loops  involving  the  new  node  Reduced  costs  of 
the  other  arcs  adjacent  to  the  new  node  remain  the  same  as  immediately  before  contract 
Note  that  we  have  at  most  one  contrax:ted  node  at  any  point  during  eliminate-chain,  but 

contracted  nodes  can  be  nested. 

The  UNCONTRACT(.s.)  operation,  applied  to  a  contracted  node  s,  restores  the  graph  as  it 
was  just  before  the  corresponding  CONTRACT  operation  and  adds  p{v)  to  prices  of  all  no  es 
in  Si.  At  the  end  of  the  chain  elimination  process,  we  apply  uncontract  until  the  oiiginal 

graph  is  restored. 

Contraction  is  used  for  efficiency  only  and  does  not  change  the  price  function  computed  by 
ELIMINATE-CHAIN,  because  by  Lemma  6.1  S,  C  Sj  for  1  <  i  <  J  <  <• 

Second  we  implement  the  search  for  the  nodes  reachable  from  Wi's  in  the  admissible  graph 
i„  Awaysiinila,  Dial’s  in,ple,nentatio„  [4]  of  Dijkstra’s  algorithm. 

a  nrioritv  queue  that  holds  items  with  integer  key  values  in  the  range  [0, . . . ,  37i],  the  amoitized 
^osl  oole  pticity  <,ueus  operations  is  constant.  We  assume  the  following  <,ueue  operauons. 

•  enqueue{v,Q):  add  a  node  u  to  a  priority  queue  Q. 

•  min{Q):  return  the  minimum  key  value  of  elements  on  Q. 

•  extract-min{Q):  remove  a  node  with  the  minimum  key  value  from  Q. 

•  decrease-keij{v,x):  decrease  the  value  of  key(u)  to  .r. 

•  .shifliQ.S):  add  (S  to  the  key  values  of  all  elements  of  Q. 

All  of  these  operations  except  Mft  are  standard;  a  constant  tinre  implementation  of  Mt  is 
trivial. 


1  procedure  scan(t;); 
for  all  (v,tjy)  do 

if  key(w)  =  oo  then 
mark  w  as  labeled; 
k€y{w)  lp{v,w); 
insert{w,Q); 

else  if  w  is  labeled  and  key{w)  <  h(lp{v,w))  then 
decr€ase-k€y{w  y  lp{v,w))] 
mark  v  as  scanned; 

I  end.  _ _ _ _ _ _ 


Figure  2:  The  scan  operation. 


Note  that  if  p  is  (2e)-feasible  and  l,{a)  >  2ne,  then  a  can  be  deleted  from  the  graph.  We 
assume  that  such  arcs  are  deleted  as  soon  as  their  reduced  costs  become  large  enough. 

We  define  the  key  assignment  function  h  that  maps  reduced  costs  into  integers  as  follows. 


h{x)  = 


0  if  X  <  0 
[1]  otherwise. 


During  the  chain  elimination  computation,  each  node  is  unlabeled,  labeled,  ov  scanned. 
UnlS^eled  nodes  have  infinite  keys;  other  nodes  have  finite  keys.  The  pr.onty  querie  Q  conta  n 
T  11  11  arp  nnlabeled  At  the  beginning  of  iteration  t,  key{wt)  is  set 

labeled  nodes.  In, ..alkali  ke,  value  of  the  queue 

executed.  This  concludes  iteration  i. 

Next  we  prove  correctness  of  the  implementation. 

Lemma  6.2  The  sets  Si  are  computed  correctly  for  every  i  =  1, . . . ,  t- 

Proof.  For  convenience  we  define  So  =  «■  Consider  an  iteration  i.  It  is  enough  show  that  Si 

is  correct  M  l<i  <t  and  Si-i  is  correct.  ,  , ,  r  •  tUo 

,  I  rTo  n  with  the  zero  kev  value.  We  claim  that  v  is  reachable  from  u);  m  the 

Let  i>  e  a  no  consider  two  cases.  If  v  was  a  node  on  Q  with  zero  key 

current  admissible  graph,  lo  seetn  ,  cKoKIo  finm  iw  bv  Lemma  6.1.  Otherwise, 

value  at  the  beginning  of  ^  TL^canned.  We  can  make  an  inductive  assumption 

thTL"abrfr«  r  r; of  h,  M-.)  =  0  i-PH-  ^ 

therefore  v  is  reachable  from  'U/',.  i  u  »• 

ut  r  be  an  admissible  path  originating  a.  It  is  el^.v  to  see  by  induction  on  the  number 

of  arcs  on  V  that  all  nodes  on  T  are  scanned  and  added  to  b.. 


It  follows  that  at  the  end  of  iteration  i,  5,  contains  all  nodes  reachable  from  Wi  in  the 
admissible  graph.  ■ 

Lemma  6.3  eliminate-chain  runs  in  0{m)  time. 

Proof  Each  node  is  scanned  at  most  once  because  a  scanned  node  is  marked  as  such  and 
Tevefadld  to  Q.  A  contracted  node  is  never  scanned.  The  t.me  to  scan  a  (noncontracted) 
node  is  proportional  to  degree  of  the  node,  so  the  total  scan  time  is  (m). 

..  •  /^(t  j.  4-  m'\  where  n^  is  the  number  oi  nodes 

The  time  of  a  contract  arcs  between  these  nodes.  The  number  of  contract 

being  contrMle  an  m  is  «  "  contract  operations  is  at  most  m  and  the 

:r“  ie:  ”■  W::Tn  T^s  the  tot.  cost  of  contract  operations  is  0(m). 

*ic  where  sind  a.re  the  same  as 

in  n:“crora«o~.  This  the  totai  time  for  the.  operations  is  0,m,, 


7  Faster  Algorithm 


•  f  oarinrp  an  Of  /umlos  N)  algorithm  for  finding  a  feasible  price  function. 
In  this  section  we  intioduce  an  0(v/  g-  )  , reduces  k  by  at  least  ^/k  at  each 

r ire:r  Sri  res..  „f  ectio„s  ^  C  ai. 

the  foliovving  iemma,  ivhich  is  related  to  Dilworth’s  Theorem  (see  e.g.  (h]). 

Lemim  7.1  Suppose  G,  <s  Uiat^Ming  S  reducer  the  number  of 

i:  ^'“'/.mtwo-.,  snch  S  con  he  found  in  0(m)  time. 


Proof.  Define  a  length  function  /'  on  Ep  by 


/'(«)  - 1  0  ( 


-1  if  a  is  improvable 


otherwise. 


The  absolute  value  of  the  path  length  »ith  respect  to  1'  is  equal  to  the  number  of  improvable 
aics  on  the  path.  ^ 

Add  a  source  node  i  to  and  .  y  R  .ive  the  shortest  paths  distances 

;romt»m'tp«:o  i  ll:  G'  sLTg-  .  acycnc'd'  can  be“computed  in  linear  time.  Dehne 

°  if' D  >  /r.  then  a  shortest  path  from  r  lo  a  node  o  with  d'(f)  =  -O  contains  a  chain  with 

at  least  \/k  improvable  arcs. 


procedure  REFlNE(e,p); 

k  the  number  of  improvable  nodes; 
repeat 

DECYCLE(Gp); 

5  4—  a  chain  or  a  set  as  in  Lemma  7.1; 
if  5  is  a  chain 

1  eliminate-chain(5); 

else 

cut-relabel(5); 

k  4-  the  number  of  improvable  nodes; 
until  A:  =  0; 
return(p); 

end. . . 


Figure  3:  An  efficient  implementation  of  REFINE. 


If  D  <  \/fc,  then  the  partitioning  of  the  set  of  improvable  nodes  according  to  the  value  of 
d'  on  these  nodes  contains  at  most  Vk  nonempty  subsets.  Let  X  be  a  subset  containing  the 
maximum  number  of  improvable  nodes  and  let  i  be  the  value  of  d’  on  X.  Observe  that  A 
contains  at  least  \/k  improvable  nodes.  Define  5  =  {u  €  V\d'{v)  <  i}. 

Clearly  Y  C  5.  Also,  5  is  closed.  This  is  because  if  «  €  5  and  there  is  a  path  from  ti  to  u; 
in  Gp,  then  the  length  of  this  path  with  respect  to  V  is  nonpositive,  so  d'{w)  <  d'{v)  <  i  and 

therefore  to  £  S. 

We  show  that  after  cut-RELABEL  is  applied  to  5,  nodes  in  X  are  no  longer  improvable.  Let 
.T  €  -Y  and  let  (u.  x)  be  an  improvable  arc.  Then  l'{v,  x)  =  -1  and  therefore  d  {v)  >  d  (x)  =  i. 
Thus  V  ^  S  and  [v,  ir)  is  not  improvable  after  relabeling  of  S.  ■ 


The  0(v/^m)  implementation  of  REPINE  is  described  in  Figure  3.  The  implementation 
reduces  the  number  of  improvable  nodes  k  by  at  least  Vk  at  each  iteration  by  eliininating 
cycles  in  Gp,  finding  5  as  in  Lemma  7.1,  and  eliminating  at  least  Vk  improvable  nodes  in  5 
using  techniques  of  sections  4,  5,  and  6, 

Lemma  7.2  The  implementation  o/ REFINE  described  in  this  section  runs  tn  0(\/nm)  time. 


Proof.  Each  iteration  of  REFINE  take  0(m)  time  by  the  results  of  the 
iteration  reduces  k  by  at  least  Vk,  and  0{Vk)  iterations  reduce  k  by 
The  total  number  of  iterations  is  bounded  by 


previous  sections.  Each 
at  least  a  factor  of  two. 


E 


o(V^). 


Corollary  3.2  and  Lemma  7.2  imply  the  following  result. 


Theorem  7.3  The  shortest  paths  algorithm  with  REFINE  implemented  as  described  in  this 
section  runs  in  0{y/nm\ogN)  time. 


8  Tighten  Operation 

In  this  section  we  describe  an  alternative  to  the  cut-RELABEL  operation  which  we  call  tight¬ 
en  This  operation  is  motivated  by  the  operation  described  in  [13]  in  the  context  of  mm.mum 
cost  flows.  Let  p  be  a  (2€)-feasible  price  function.  Assume  that  we  eliminated  cycles  m  Gp, 
and  let  d\  and  P  be  as  in  the  proof  of  Lemma  7.1.  Define  a  new  price  function  p  by 

d'iv) 

p'{v)  =  p{v)  +  e 

TIGHTEN  computes  d!  and  replaces  p^y  P  •  This  takes  0(m)  time. 

For  any  u,  0  <  p{v)  -  p'(v)  <  6.  Thus  if  l^v,  w)  >  0  then  lA^M  >  If -^  <  | 

0  then  l'{v,w)  =  0  and  thus  d'{w)  <  d'{vy,  therefore  /p-(i;,«i)  >  />,«;)  >  £.  Finally  if 

C(v,xu)  <  -e,  then  l'{v,w)  =  -1  and  thusd'H  <  d'{v)-y  therefore  >  l,{v,w)-elD. 

This  implies  that  TIGHTEN  creates  no  improvable  arcs.  Since  D<n-h  the  reduced  cost 
of  every  existing  improvable  arc  increases  by  at  least  Therefore  the  implementation  of 
REFINE  bcised  on  TIGHTEN  takes  O(Tim)  time. 

Note  that  TIGHTEN  does  not  maintain  integrality  of  p,  so  the  method  cannot  terminate 
whe^r  rethes  1  However,  if  .  =  0(l/n)  and  /  is  integral,  an  c-feasible  price  function  can 
be  converted  into  a  feasible  price  function  using  rounding  and  a 
computation.  Therefore  the  overall  running  time  of  the  algorithm  is  0(nmlog(ni  )). 

This  bound  can  be  improved  by  using  tighten  in  combination  with  eliminate-chain. 
Implemented  this  wav,  REFINE  works  as  follows.  It  starts  by  removing  eye  es  rom  Tp 
computino-  d'  D  and  I  If  D  >  v^,  then  eliminate-chain  is  applied  to  the  appropriate  chain, 
eiimhiatiirg  at  least  v/t  improvable  nodes.  If  D  <  ^,  then  TIGHTEN  'is  applied,  increasing  the 
reduced  cost  of  every  improvable  arc  by  at  least 

The  first  case  cannot  occur  more  than  0(s/^  times  since  all  improvable  nodes  will  be 
eliminated.  The  second  case  cannot  occur  more  than  0{s/n)  times  since  reduce  cos  o  eve\ 
tprovtwe  arc  will  increase  by  at  least  e  and  these  arcs  will  no  longer  be  improvable.  The 

resulting  algorithm  runs  in  Oi^/nn^\og(nN))  time. 


9  Concluding  Remarks 

We  described  a  framework  for  designing  sraliiig  algorithms  and  Wo 

and  TIGHTEN,  that  can  be  used  to  design  algorithms  within  this  framework  Tli« 

is  very  Ilexibie  and  can  be  used  to  design  numerous  algorithms  for  the  problem.  g 

res  1  e  tp  oved  the  time  bound  for  the  problem.  We  believe  that  further  ."'  es  r^ga  m  o 

S  framework  is  a  promising  research  direction.  Our  algorithms  are  easy  to  implement  and 


9 


may  have  practical  implications;  this  work  was  in  fact  motivated  by  an  experimental  study  of 
minimum-cost  flow  algorithms  [12]. 

If  one  is  interested  only  in  the  algorithms  based  on  cut-relabel,  bit  scaling  can  be  used 
instead  of  e-feasibility.  This  version  of  the  algorithm  rounds  lengths  up  to  a  certain  precision, 
initially  the  smallest  power  of  two  that  is  greater  or  equal  to  N.  Each  iteration  of  the  algorithm 
starts  with  a  price  function  that  is  feasible  with  respect  to  the  current  (rounded)  lengths. 
At  the  beginning  of  an  iteration,  the  lengths  and  prices  are  multiplied  by  two,  and  one  is 
subtracted  from  the  arc  lengths  as  appropriate  to  obtain  the  higher  precision  the  lengths.  Te 
resulting  price  function  is  1-feasible  with  respect  to  the  current  length  function;  the  feasibility 
is  restored  using  the  cut-relabel  operations  as  described  above.  The  bounds  obtained  for  the 
algorithms  that  use  the  relabel  operation  remain  valid,  but  some  aspects  of  the  implementation 
and  analysis  become  slightly  simpler.  In  particular,  the  only  possible  negative  reduced  cost 
value  is  -1  and  the  length  function  /'  defined  in  the  proof  of  Lemma  7.1  is  the  same  as  Ip. 
Also,  no  rounding  is  required  for  computing  the  priority  queue  keys.  However,  the  tighten 
operation  cannot  be  defined  in  a  natural  way  if  bit  scaling  is  used. 

Our  definition  of  ^-feasibility  corresponds  to  that  of  ^-optimality  for  minimum  cost  flows 
[11  14].  If  one  follows  [11, 14]  faithfully,  however,  one  would  define  e-feasibility  using /p(a)  ^  f 
instead  of  (2)  and  not  consider  arcs  with  zero  reduced  costs  admissible.  Under  these  definitions, 
the  admissible  graph  cannot  have  zero  length  cycles,  so  there  is  no  need  for  decycle.  However, 
these  definitions  seem  to  lead  to  an  0(log(niV))  bound  on  the  number  of  iterations  of  the  outer 

loop  of  the  method. 

In  conclusion  we  would  like  to  mention  a  natural  variation  of  the  TIGHTEN  operation  related 
to  continuous  optimization  techniques.  Suppose  we  use  Ip  instead  of  /'  and  redefine  p  y 

p'{v)  =  p+Sd'iv). 

We  can  interpret  d'  as  the  direction  we  want  to  move  in,  and  d'  as  a  parameter  that  determines 
the  step  size.  Then  we  can  define  a  penalty  function  $  whose  value  is  determined  by  the 
reduced  costs,  and  pick  6  to  achieve  a  large  decrease  in  For  example,  vve  can  define 
absolute  value  of  the  most  negative  reduced  cost  and  set  d  —  1  -f  ■  len  an  app  ica 
TIGHTEN  reduces  $  to  at  most 

D  , 

^  0  +  *^“^  n  +  l' 

(The  last  inequality  follows  from  the  fact  that  D  <  n$.)  The  resulting  algorithm  runs  in 
0(nmlog(n[V))  time.  Using  a  different  penalty  function  may  give  a  different  result. 
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